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NLO corrections to AF = 2 effective Hamiltonians and SUSY 
contribution to FCNC 

Ignazio Scimcmi a 

a Dip. di Fisica, Universita degli Studi di Roma "La Sapienza", P.le A. Moro 2, 00185 Roma, Italy. 

The most general QCD-NLO Anomalous-Dimension Matrix of all four-fermion dimension six AF = 2 operators 
is presented. Two applications of this Anomalous-Dimension Matrix to the study of SUSY contribution to K — K 
mixing are also discussed. 



1. Introduction 

Theoretical predictions of several measurable 
quantities, which are relevant in K-, D- and B- 
meson phenomenology, depend on the matrix el- 
ements of some AF = 2 four-fermion operators. 
Examples are given by FCNC effects in SUSY ex- 
tensions of the Standard Model @-(2), by the B°- 
B® width difference (|] and by the 0(1 /ml) cor- 
rections in inclusive b-hadron decay rates (which 
actually depend on the matrix elements of several 
four-fermion AF = operators In all these 
cases, the relevant operators have the form 



Q = C aPp °{b a Tqp)(b p T q<J ), 



(1) 



where T is a generic Dirac matrix acting on (im- 
plicit) spinor indices; a-a are colour indices and 

(jocfipo ig either gafi gpa or $ a <J 5 p[3^ AU the oper _ 

ators discussed in this paper appear in some "ef- 
fective" theory, obtained by using the Operator 
Product Expansion (OPE). As a consequence, in 
all cases, three steps are necessary for obtaining 
physical amplitudes from their matrix elements: 

i) matching of the original theory to the effec- 
tive one at some large energy scale; 

ii) renormalization-group evolution from the 
large energy scale to a low scale suitable for the 
calculation of the hadronic matrix elements (typ- 
ically 1-5 GeV); 

iii) non-perturbative calculation of the hadron- 
ic matrix elements. 

In the following a calculation of the two-loop 
Anomalous Dimension Matrix (ADM) relevant 
for AF = 2 transition amplitudes and its appli- 
cation to FCNC in SUSY extension of the Stan- 



dard Model (SM) is presented. I refer to the 
original papers for all details of the calculations, 

1. e., respectively to ref. Jq| for the calculation of 
the ADM and to ref. for its applications to 
FCNC effects in Supersymmetry. 

The paper is organized as follows. In sec. || I 
introduce the operator basis, the Wilson coeffi- 
cients and the Renormalization Group Equations 
(RGE). In sec. | constraints on SUSY parame- 
ters coming from FCNC in the K — K mixing are 
derived. 

2. Effective Hamiltonian formalism and 
ADM 

The three steps (i-ii-iii, sec.|l|) needed to use 
OPE are treated in the following. 

In all cases of interest, the matrix elements of 
the effective Hamiltonian can be written as 



(F\H eff \I) = J2( F \Qi(»)\ I )Ci(»): 



(2) 



where the Qi(^)s are the relevant operators renor- 
malized at the scale /x and the Ci(^t)s are the cor- 
responding Wilson coefficients. 

The set of operators relevant for a description 
of FCNC effects in SUSY models is: 



Q 2 

Q:>, 

Qa 

0.5 



d a ^(l- l5 ) S a J% t (l- 75 )A 

d a (i- l5 ) s a <F(i- 75 )A 

J«(l_ 75 ) s /3 J/3(l_ 75 ) s «, 

d a (l- l5 )s a ^(1 + 75 )A 
<?*(!- TsV 3 dP (1 + 75)5°, 



(3) 



together with operators <3i,2.3 which can be ob- 
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tained from Qi 2 3 by the exchange (1 — 75) <-> 
(1+75). 

Let us now represent the operators as row vec- 
tors Q and the coefficients, C{n), as column 
ones. The vectors C{p) are expressed in terms 
of their counter-part, computed at a large scale 
M, through the renormalization-group evolution 
matrix W[fi, M] 



C{n) = W[ii,M]C{M) 



(4) 



The initial conditions for the evolution equations, 
C(M), are obtained by matching the full theory, 
which includes propagating heavy-vector bosons 
(W and Z°), the top quark, SUSY particles, etc., 
to the effective theory where the W, Z°, the top 
quark and all the heavy particles have been re- 
moved. In general, C(M) depend on the defini- 
tion of the operators in a given renormalization 
scheme. 

The coefficients C(/z) obey the renormaliz- 
ation-group equations: 



l T { a s) 



C(t,a s (t),X(t)) = 0, (5) 



where t = ln(M 2 //i 2 ). The term proportional 
to P\, the (5- function of the gauge parameter 
\{t) (for covariant gauges), takes into account the 
gauge dependence of the Wilson coefficients in 
gauge-dependent renormalization schemes, such 
as the RI scheme |n|. This term is absent in 
standard MS schemes, independently of the regu- 
larization which is adopted (NDR, HV or DRED 
for example). 

If no quark threshold between the scales M and 
/i is crossed, eq. (^) can be easily solved (the 
formulae for the general case can be found in 
refs. 11 12j]). At the next-to- leading order, one 
can write 

W\p, M] = M[fj]U\p, M]M -1 [Af] , (6) 
where U is the leading-order evolution matrix 



U\p, M] 



a.{M) 



/2/3o 



(7) 



and the NLO matrix is given by 
« s (m) 



1 



4-7T 



- J[A( M )] • 



(8) 



By substituting the expression of the C(n) 
given in eq. (ji]) in the renormalization-group 
equations (|J), and using W\p, M] written as in 
eqs. (||)-(||), one finds that the matrix J satisfies 
the equation 



J- 



01 x 



dJ 
9A 



J, 



7 



(0)T 



2/?o 



,(0)T. 



7 



(1)T 



2/3o 



.(9) 



In eqs. and (||), /3 , f3\ and /?" are the first 
coefficients of the /^-functions of a s and of A, re- 
spectively; 7^°^ and are the LO and NLO 
anomalous dimension matrices. U is determined 
by the LO anomalous dimension matrix 7^) and 
is therefore regularization and renormalization- 
scheme independent; the two-loop anomalous di- 
mension matrix 7*- 1 -* , and consequently J and 
W\p, M], are, instead, renormalization-scheme 
dependent. The precise definition and a complete 
discussion of the scheme dependence of eq. (j^) can 
be found in ref. || and in references therein. 

Eventually, the B-parameters which paramet- 
rize the hadronic matrix elements of eq. (^) have 
been evaluated non-perturbatevely on the lat- 
tice HfuJ. 

In what follows I have used the extimate of 
ref. @. 

3. FCNC in SUSY 

The considered operator basis is the one rele- 
vant for the analysis of FCNC in SUSY in the so- 
called mass insertion approximation ||. In this 
approach for fermions and sfermions states all the 
couplings of these particles to neutral gauginos 
are flavor diagonal and FCNC effects are shown 
by the non-diagonality of sfermion propagators. 
The pattern of flavor change, for the K-system, is 
given by the ratio 



(4) 



(mi) 2 



AB 



ijJAB 



M 2 



(10) 



where {rnf^)\ B are the off-diagonal elements 
of the d mass squared matrix that mixes fla- 
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vor i, j for both left- and right-handed scal- 
ars (A, B =Left, Right \ and M sq is an average 
squark mass, see e.g. Q. The sfermion propa- 
gators are expanded as a series in terms of the 
5's and the contribution of the first term of this 
expansion is considered. Assuming that all s- 
particles are heavier than fermions one produces 
effective four fermions interactions and OPE can 
be used. 

The Wilson coefficients at the matching scale 
M sq can be found, e.g., in ref. and I refer 

to these articles for their complete expression. 

One can provide a set of constraints on SUSY 
variables coming from the Kl — K$ mass differ- 
ence, AMk and the CP violating parameter ck 
defined as 



AM K = 2Re(K°\H eS \K°), 
1 



(K 



V2AM K 



lm(K \H cS \K ) 



(11) 



The parameter space is composed of two real and 
four complex entries, that is M sq , rrig and (Si 2 )ll, 

(«?a)HR. 

Neglecting interference among different SUSY 
contribution one gives upper bound on the (5's, 
giving specific values to M sq and rrig. I consider 
the conditions M, 



S(l 



rrig in subsec. 



3.1 



and the 

case in which the scalars of the first-two genera- 
tions are heavier than the rest of the supersym- 
metric spectrum in subsec. 3.2. Constraints on 



individual (5's are provided. Indeed, it is meaning- 
ful to study the interference of cancellation effects 
only in specific models. 

The physical condition used to get the bounds 
on the 5' s is that the SUSY contribution (propor- 
tional to each single S) plus the SM contribution 
to AMk an d e K do not exceed the experimental 
value of these quantities. 



sq 



3.1. M t 

For Left-Right mass insertions, I consider two 
possible (extreme) cases: \{8f 2 )LR\ ^> \(Sf 2 )RL\ 
and \{5f 2 )LR\ = \[$v2)rl\- I n the second case the 
contributions proportional to {8f 2 ) 2 LR , {Si 2 ) RL 
and (Si 2 )lr(5i 2 )lr are combined. This approach 
improves the one of refs. p,^5[, where contribu- 
tions proportional to {8f 2 )LR.{5f 2 )LR were consid- 
ered independently from the ones proportional to 



(6f 2 ) 2 LR and (Sf 2 ) RL . 

(Si 2 )ll(Si 2 )rr contribution can be treated in- 
dependently from the {8f 2 ) 2 LL and {Sf 2 ) 2 RR ones, 
since the latter ones do not generate Left-Right 
operators and they are suppressed with respect 
to the first contribution. 

Results for the bounds of the S's for M sq ~ 200 
GeV and various values of the ratio x — to? /M 2 q 
derived for real parts of the 5's are shown in 
tab. 0. Let us compare these results with the 
previous analyses of refs. In both papers 

the matrix elements were computed in the Vac- 
uum Insertion Approximation (VIA) and with no 
evolution [Q or leading order evolution || of the 
coefficient functions. In this analysis lattice B- 
parameters and NLO evolution of the coefficients 
have been calculated. This allows a consistent 
matching of the renormalization-scheme and /z 
dependence of the QCD-corrected TtffjT 2 . 

3.2. Heavy first-two generation scalars 

Let us assume now that the scalars of the first 
two generations (the typical mass of which is here 
again denoted by M sq ) is much heavier than the 
rest of the SUSY spectrum (which has a typi- 
cal mass of rrig). In this case the heavy squarks 
have to be integrated out as a first, at a scale 
M, 



sq- 



Then, the coefficients must be evolved to 
the scale m g where all the rest of the supersym- 
metric particles are also integrated out. 

In table || results are shown for the real parts 
of the 5's with nig ~ 250 GeV and several values 
of M sq . 

The constraints that come from the four pos- 
sible insertions of the S's are presented: in the 
first and second rows only terms proportional re- 
spectively to {5i 2 )ll and (Sf 2 )LR are considered; 
in the last two rows the contribution of opera- 
tors with opposite chirality, RR and RL, is also 
evaluated by assuming (Si 2 )ll — ($12) RR an( i 
(^i 2 )l_r = (Su)rl- 

The theoretical improvements of this results 
with respect to ref. |J have already been discussed 
in subsec. 3.1. The combination of B-parameters 



and NLO-QCD corrections change the LO-VIA 
results by about 25 — 35%. As expected Q, the 
tightest constraints are for the cases (<5f 2 )iL = 
{Sf 2 )RR and (Sf 2 ) LR = (6( 2 )rl- In these cases 
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No-QCD, 


LO, 


LO, 


NLO, 




VIA 


VIA 


Bi 


Bi 




xl0~ 3 


xlO" 3 


xlO" 3 


xlO" 3 


X 




0.3 


5.0 


5.7 


7.7 


7.7 


1.0 


11 


12 


16 


16 


4.0 


25 


29 


39 


39 


X 


VIM^lrI {\(S&)lk\ » \(St 2 ) RL \) 


0.3 


1.1 


0.84 


1.1 


0.96 


1.0 


1.2 


0.93 


1.2 


1.1 


4.0 


1.8 


1.3 


1.6 


1.5 


X 


v/|Re(^ 2 




= {$12) Kl) 


0.3 


2.0 


1.4 


0.89 


0.67 


1.0 


1.1 


0.97 


1.8 


3.0 


4.0 


1.3 


1.0 


1.4 


1.3 


X 




ll(<5i 2 )rr| 


0.3 


0.64 


0.39 


0.40 


0.33 


1.0 


0.71 


0.44 


0.45 


0.37 


4.0 


1.0 


0.61 


0.62 


0.52 



Table 1 

Limits on Ke(Sij) AB {Sij) CLn with A,B,C,D = 
(L,R), for an average squark mass M sq = 200 
GeV and for different values of x = rri~ /M s 2 ? . 



M aq 


No-QCD, 


LO, 


LO, 


NLO, 


[TcV] 


VIA 


VIA 


Bi 


Bi 




xlO" 2 


xlO -2 


xlO" 2 


xlO -2 






2 


3.1 


3.6 


4.9 


4.9 


5 


7.5 


8.8 


0.12 


0.12 


10 


15 


18 


25 


24 






Re(5f 2 )| K | 




2 


2.1 


1.4 


1.8 


1.6 


5 


9.8 


6.5 


8.2 


7.2 


10 


34 


22 


28 


25 




V|Re(S? 2 


) LR = Re(Sf 2 ) RL \ 




2 


0.66 


0.35 


0.35 


0.28 


5 


1.5 


0.77 


0.79 


0.64 


10 


3.0 


1.5 


1.5 


1.2 




V]Re(5f 2 


)L=Re(««l 




2 


1.1 


0.52 


0.51 


0.41 


5 


4.1 


1.6 


1.6 


1.3 


10 


10 


3.6 


3.4 


2.7 



Table 2 

Limits on Ke(Sf 2 )AB from AMk with gaugino 
masses of 250 GeV. 



the coefficients proportional to (Sf 2 ) ll {5f 2 ) RR-, 
12) l r(S f 2) rl dominate the others. 
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